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ABSTRACT 

Recent radial-velocity surveys for GK clump giants have revealed that planets also exist 
around ~1.5— 3 M© stars. However, no planets have been found inside 0.6 AU around clump 
giants, in contrast to solar-type main-sequence stars, many of which harbor short-period planets 
such as hot Jupiters. In this study we examine the possibility that planets were engulfed by host 
stars evolving on the red-giant branch (RGB). We integrate the orbital evolution of planets in 
the RGB and helium burning (HeB) phases of host stars, including the effects of stellar tide and 
stellar mass loss. Then we derive the critical semimajor axis (or the survival limit) inside which 
planets are eventually engulfed by their host stars after tidal decay of their orbits. Especially, we 
investigate the impact of stellar mass and other stellar parameters on the survival limit in more 
detail than previous studies. In addition, we make detailed comparison with measured semimajor 
axes of planets detected so far, which no previous study did. We find that the critical semimajor 
axis is quite sensitive to stellar mass in the range between 1.7 and 2.1 Mq, which suggests a need 
for careful comparison between theoretical and observational limits of existence of planets. Our 
comparison demonstrates that all those planets are beyond the survival limit, which is consistent 
with the planet-engulfment hypothesis. However, on the high-mass side (> 2.1M Q ), the detected 
planets are orbiting significantly far from the survival limit, which suggests that engulfment by 
host stars may not be the main reason for the observed lack of short-period giant planets. To 
confirm our conclusion, the detection of more planets around clump giants, especially with masses 
> 2.5Mq, is required. 

Subject headings: planetary systems: evolution — stars: evolution 



1. INTRODUCTION 

Detection of a significant number of exoplan- 
ets orbiting Sun-like stars (FGK dwarfs) has 
helped us to improve our understanding of the 
origins and diversity of planetary systems. Statis- 
tics on detected planets suggests the presence 
of several different populations and correlations. 
That has motivated theorists to synthesize planet 
populations by integrating the formation process 



from planetesimals to nlanets (llda & Lin 


2004a 


1). 


2005.l2008allb. 201C* Kennedv & Kenvon 


2008a 


1): 


Mordasini et al.l 2009allb). Cormoarison with sta- 
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tistical properties from observations has then en- 
abled us to calibrate theories for planet formation. 

As for heavier stars such as BA dwarfs, the 
situation is somewhat different. Precise radial- 
velocity measurement is not successful, because of 
the lack in absorption lines on their hot photo- 
spheric surfaces. Instead, efforts have been made 
to find planets around GK giants that used to be 
BA dwarfs on their main-sequence. Thus, com- 
parison between observations and theories is not 
straightforward, unlike in the case of FGK dwarfs. 
That means we have to take into account modifi- 
cations to orbital configurations of planets during 
their host stars' evolution. 

According to recent radial-velocity surveys for 



GK clump giants, it appears th at there is a lack 



of giant planets inside 0.6 AU ( Sato et al. 2010L 



and references therein). Such a deficit is not seen 
around Sun-like stars which often harbor short- 
period planets such as hot Jupiters. It is un- 
clear whether the deficit is primordial — namely, 
short-period giant planets (SPGPs) are not formed 
around BA dwarfs originally; otherwise short- 
period planets are removed during host stars' evo- 
lution. The latter possibility is focused on in this 
paper. 

A star evolves off the main sequence (MS) to- 
wards the red-giant branch (RGB), after exhaus- 
tion of hydrogen at its center. Then, once the 
helium ignites, the star enters a next central nu- 
clear burning phase that is sometimes called the 
helium-burning phase (HcB). Most of the clump 
giants with detected planets are thought to be on 
their HeB. Before reaching the HeB, namely, in 
the RGB phase, the star expands substantially, 
becomes highly luminous, and may lose substan- 
tial mass, which should affect orbits of surround- 
ing planets. 

The evolution and fate of planets in the RGB 
phase of their host stars have been investigated in 
the context of the surv ival of the E arth. Based on 
his simple estimation, Vilal ( 19841). followed by a 
somewhat detailed argument bv lGoldstein ( 1987 ). 
showed that the Earth would lose its momentum 
quickly in the extending envelope of the Sun before 
the envelope starts to shrink again. That means 
the Earth would be unable to escape from the Sun 
once engulfed. Those studies suggest the need for 
the detailed simulation of the Sun's evolution. The 
Sun's evolution, incl uding solar mass l oss, wa s sim- 
ulated in detail by ISackmann et al. (1993), who 
demonstrated that the decrease in the solar mass 
is effective in pushing the terrestrial planets other 
than Mercury outwards and preventing them from 
being engulfed by the expanded Sun. Later, build- 
ing on theories for formation of close b inaries (e.g. 



Livio 



1982ULivio fc Sokerlll983lll984l) . lRasio et al 



(|1996I ) incorporated the tide raised on the Sun 
and then demonstrated that Venus is likely to fall 
into the Sun because of the solar tide enhanced by 
the expansion of the envelope, while the Earth's 
fate remained inconclusive, depending on uncon- 
strained parameters for tide. Other effects such 
as the gr avitational drag, ram pressure, wind- 



and orbital instability due to secular perturbation 
(jDuncan fc Lissauerlll998l ) are known to be negli- 
gible. Note that the last possibility could be im- 
portant in extrasolar planetary systems, depend- 
ing on orbital configurations. 

To understand the obser ved lack of SPGPs 
arou nd GK clump giants , ISato et al.l (|2008aj ) 



and IVillaver fc Livid ([iobH ) investigated the or- 
bital evolution of planets around i ntermediate- 



mass RGB stars in a similar way to iRasio et al 
(|l996l ). They demonstrated that SPGPs did un- 



dergo orbital decay due to stellar tide, ending 
up bein g swallowed by t heir h ost stars. In par- 
ticular, Villaver fc Liviol ( 20091 ) derived a critical 
scmimajor axis, a cr it (note that they called it a 
minimum scmimajor axis and denoted it by a m i n ), 
beyond which planets evade engulfmcnt by their 
host stars during the RGB phase: In four cases 
of different stellar initial masses, M^^ (i.e., 1, 2, 
3, and 5 M Q ), the values of a crit that they de- 
rived for Jovian-mass planets were 3.0, 2.1, 0.18, 
and 0.45 AU, respectively. Then, without detailed 
comparison of a cr it with the orbital distribution 
of known planets around GK clump giants, they 
concluded that engulfmcnt by host stars may be 
responsible for the observed lack of SPGPs inside 
0.6 AU and any ad hoc mechanism to never make 
SPGPs may not be required. 

In this paper, to verify the planet-engulfment 
hypothesis, we make a detailed comparison of 
the theoretical limit of survival (namely, a cr it) 
with semimajor axes of the planets detected 
so far around GK clump giants. The mass of 
those giants ranges from ~ 1.5 and 3 M@. As 
known in the study of stel lar evolution (e.g., 
Kippenhahn fc WeigertJ Il99dh . the RGB-tip ra- 



dius is sensitive to stellar mass in such a range. 
Given that stellar tide depends strongly on stel- 
lar radius (see section [2]), the limit of survival 
from engulfment by host stars should be sensi- 
tive to host stars' mass, which implies the need 
for detailed comparison between the theoretical 
and observational limi ts. It should be noted that 
Nordhaus et~aTl (|2010h has recently made a thor- 



accretion (Fviial Il984t iDuncan fc Lissauerl LL998), 



ough study of the tidal evolution of orbits of plan- 
ets around post-MS stars. Their focus was, how- 
ever, on planets around white dwarfs that survive 
engulfment by AGB stars. Of interest in this pa- 
per is planets around RGB or HeB stars that have 
not entered the AGB phase yet. 
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This paper is organized as follows. In section [2] 
we first describe our physical model and computa- 
tion method to simulate the orbital evolution of a 
planet around its host star that is evolving off the 
MS. In section [3J we derive a critical semimajor 
axis beyond which planets survive the host stars' 
RGB /HeB phase, and investigate its sensitivity to 
stellar mass and other parameters. In section HI 
we compare the derived survival limit with mea- 
sured semimajor axes of the detected planets. The 
comparison suggests that planet engulfmcnt may 
not be the main reason for the observed paucity 
of SPGPs, apart from statistical sufficiency. In 
section [5j we evaluate the impacts of several un- 
certainties on the survival limit and find that our 
conclusion is not affected by them. In section|6l we 
mention the possibility of stellar-mass-dependcnt 
formation process that may account for the ob- 
served paucity of SPGPs. Finally, we summarize 
this paper, claiming the need for more planet sam- 
ples around giants of > 2.5M Q to confirm the va- 
lidity of our findings, in section [7] 

2. PHYSICAL MODEL AND COMPU- 
TATION METHOD 

We simulate the evolution of orbits of planets 
during host stars' evolution. We assume circular 
orbits, because our focus is on semimajor axis in 
this study. The eccentricities of planets detected 
around the GK clump giants that we target in this 
study are as small as 0.25. Such a small eccentric- 
ity has a small effect on the evolution of semima- 
jor axis. The impact of the finite eccentricity on 
the limit of survival is evaluated in section I 



Note that previous studies ( Villaver fc Livioll200 



E3 

!009: 



Nordhaus et al.ll2010l ) also assumed zero eccentrit- 
icies. 

We consider the effects of stellar tide and 
mass loss on the planetary orbit. We neglect 
other competing processes such as the frictional 
and gravitational drag forces by stellar wind and 
change in the planet mass due to stellar-wind ac- 
cretion and evapor ation, which wer e eval uated 
to be negligible by Villaver fc Livid ( 20091 ) and 
Duncan fe Lissauen ( 19981) . Thus, we integrate 
the equation, 



Ida 
a dt 



= -6; 



k M p 



1 



Mp 



(1) 



a given mass. In the above equation, a is the semi- 
major axis, t is time, M p is the planet's mass, M* 
and R+ are the mass and radius of the host star, 
respectively, k is the apsidal motion constant, and 
T is the eddy turnover timescale (see below). The 
first term on the rig ht-hand si de represents the ef- 
fect of stellar tide (|Hutlll98ll ). Since RGB stars 
are slow rotators, we assume no stellar rotation, 
which means that the stellar tide always causes 
orbital decay of the planet. 

As for the paramete rs for stella r tide such 
as k and T, fo l lowin g Rasio et al. ( 19961 ) and 
Villaver fc Liviol (|2009l ). we adopt the turbulent 
viscosity: 

where M onv is the mass in the convective envelope 
and / is given by 



P 
2T 



(3) 



with P being the orbital period. The eddy 
turnover timescale T is given by 



T 







1/3 



(4) 



where R cm is the radius at the base of the convec- 
tive zone and L* is the stellar intrinsic luminosity. 
The factor / is introduced to weaken the stellar 
tide when the orbital period is shorter than the 
eddy turnover timescale. However, it is a matter 
of debate how to deal with such fast tide (e.g., 
IZahnll2008l ). We discuss the impact of the uncer- 
tainties in / on the survival limit in section 15.11 

The second term on the right-hand side of equa- 
tion (JTJ) represents orbital migration due to stellar 
mass loss (M* < 0). We use the Reimers' param - 
eterisation for stellar mass loss (|Reimersl Il975l ) , 
namely, 



M* = -4xl(T 13 ?7 



R* 



M 



©yr 

(5) 

where r\ is the mass-loss parameter of order unity. 

In this study, we simulate stellar evolution di- 
rectly with the code MESA ver. 2258 (jPaxton et all 
201lh to calculate M*(t), L*(t), M*(t), 



from the zero-age main-sequence of the host star of 



M env (t), and R cnv (t). We adopt the mass- loss 
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parameter 77 = 0.6 and the mixing-length param- 
eter a — 1.5. We include the effects of connec- 
tive overshooting, following the MESA's prescrip- 
tion where the overshoot parameters are set to 
be 0.0128 at the bottom of non-burn regions and 
0.014 at all the other boundaries. Several different 
stellar-evolution models that are simulated with 
different numerical codes are availab l e in t he lit- 
erature and online. Villaver fc Livid ( 20091 ) used 
a stellar-evolution code other than MESA. We 
discuss the differences and their impacts on the 
survival limit in section [5.21 



3. THEORETICAL LIMIT OF SUR- 
VIVAL 

3.1. Orbital Evolution 

We derive the theoretical limit of survival from 
engulf mcnt by evolving host stars. Figure [T] shows 
examples of the orbital evolution of planets around 
stars with initial masses, M*^, of (a) 1.8M Q and 
(b) 2.0Af Q ; the metallicity is 0.02. We have inte- 
grated the orbit of a lMj planet {My. Jupiter's 
mass) with different initial semimajor axes (solid 
lines) from the host-star's zero-age main sequence 
until the planet's semimajor axis falls below the 
stellar radius (a < i?*) or until the AGB ther- 
mal pulse (AGB-TP) occurs. The evolution of the 
stellar radii is represented by the dashed lines. 

Although we have continued the simulations af- 
ter the end of HeB (grey areas in the figure), our 
interest in this study is in whether planets survive 
the RGB and HeB phases of their host stars. This 
is because the planet-harboring GK clump giants 
that we target for the comparison of our survival 
limit with the semimajor axes of their planets are 
thought to be in the HeB or RGB phases, as stated 
in Introduction. A s for orbital evolution o f planets 
around AGB stars, Nordhaus et al. ( 2010h made a 
thorough investigation. 

As seen in Fig. [TJt, the 1.8M Q star expands 
up to 0.5 AU at the RGB-tip. As the stellar 
radius becomes large, tidal decay is rapidly en- 
hanced (eq. p]), so that the star consequently 
swallows planets whose initial semimajor axes, 
ttj, are smaller than 1.1 AU. The planet that 
starts out at 1.1 AU is tidally pulled by the 
host star, but it barely survives the RGB-tip 
and remains in the HeB phases. Outer plan- 
ets (e.g., a, = 2.0 AU) are found to be pushed 
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Fig. 1. — Orbital evolution of 1- Jupiter-mass planets 
with different initial semimajor axes (solid lines) for 
M* = (a) 1.8 M Q and (b) 2.0 M Q ; the metallicity, 
Z*, is 0.02 in both cases. Dashed lines represent the 
evolution of stellar radii. 



outwards because of stellar mass loss; the effect 
is, however, small, unlike in the case of IMq 



stars, as also noted previously (|Sato et al.ll2008a ; 
Villaver fc Livid l2009t iNordhaus et al.ll2010l ). As 
for the 2Mq star, the RGB-tip radius is as small 
as 0.20 AU (Fig. [lb). As a consequence, even a 
planet starting out at 0.36 AU can avoid engulf- 
ment during the RGB/HcB phases. 

We define a critical initial semimajor axis (a cr it) 
below which planets end up being engulfed by 
their host stars at some point on the RGB or HeB. 
For the above two examples, a cr ; t = 1.1 AU around 
the 1.8Af Q star, while et cr it = 0.36 AU around the 
2M Q star. It is thus illustrated that the survival 
limit is rather sensitive to the host star's mass. 
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Star's RGB-tip radius 



1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3 3.1 3.2 3.3 
Host Star's Initial Mass [solar mass] 

Fig. 2. — The critical initial semimajor axis, a cr it, 
inside which planets are engulfed by host stars before 
reaching the end of the helium-burning phase for two 
choices of planetary mass, M p , namely, M p = 1 Mj 
and 20 Mj. For comparison, the RGB-tip radius of 
the host star is shown by the dashed line. The stellar 
metallicity is 0.02 in the calculations. 



Furthermore, comparing Fig. [Tfc> with the upper - 
right panel of Fig. 1 of IVillaver fc Livid (|2009l ). 
one notices that the value of a cr it derived here for 
the 2M@ star is considerably diff erent from that 
derived bv lVillaver fc Livid (|2009f ) , who estimated 
a crit to be 2.1 AU, instead of 0.36 AU. The reason 
for the difference is discussed in section (521 



For comparison with semimajor axes of the de- 
tected planets, the above definition of a cr it may be 
inappropriate, because it is the initial value of the 
semimajor axis of the planet that barely survives 
the host-star's RGB-tip. As shown in the above 
figure, the planet's semimajor axis after the RGB- 
tip is lower than its initial value. However, the 
probability that such planets are detected should 
be low. Because of unambiguity of definition, we 
decide to use the critical semimajor axis defined 
above. 

3.2. Critical Semimajor Axes 

Figure [2] shows the relation between a cr ;t and 
M M for M p = lMj and 20 Mj. The stellar RGB- 
tip radius (i?ti P ) is also shown by the dashed line. 

First, one notices that the curves for a cr it are 
similar in shape to that for i?ti P - This is be- 
cause the tidal decay rate of planetary orbit is 



quite sensitive to the stellar radius (sec eq. [Tj). 
That means planet engulfment, in most cases, oc- 
curs right around the RGB-tip of the host star, 
at which the stellar radius takes a maximum dur- 
ing RGB/HeB phases. To see that more clearly, 
we derive an approximate expression of a„it as fol- 
lows. From equation ([TJ, a typical decay timescale 
(Tdecay) can be expressed by 



7"decay 



6k M„ 



Ik 

a 



(6) 



provided the mass-loss effect is negligible. If Tdecay 
is comparable to the typical duration of the RGB 
phase (trgb), then the planet is engulfed by its 
host star. Thus, a cr ; t would be expressed as 



^crit — 



G 



1/5 ^4/5 „ , r l/5 

J-cff M P ,1/5 R 8/5 
3/5 ^RGB^tip 



0.11 



AC 
, 8/5 /M, 



0.1AU/ \2M P 



/ TtGB \ 

\5Myv) 



1/5 



Mp 
Mj 



1/5 



-3/5 



Toff \ 

4000K / 



1/5 



AU 



(7) 



In deriving this equation, based on our numer- 
ical results, we have assumed that JW cnv ~ M* 
and i?env ~ in equations © and ^ and used 



T ~ (A4i? t 2 ip /3 J L tip ) 1 / 3 with the stellar RGB-tip 
47^7*^ and k ~ f/6 



luminosity i t ; p 

= (P/2T) 2 /6 (see eqs. [2]-[4]). Equation © con- 
firms that i?ti P is the most important factor and 
the others have minor impacts on a cr ;t. Note that 
we have found that trgb is not sensitive to stel- 
lar mass; for example, the periods during which 
> 0.5i?ti P in the RGB phase are 4.3 Myr, 
9.7 Myr, and 5.2 Myr for 1.5M Q , 2.0M Q , and 
3. OA/0 stars with = 0.02, respectively. It is 
worth noting that equation ((7|) agrees well with 
the numerical values of a cr j t on the high-mass 
side in Fig. the differences are less than 10 % 
for AI± > 2.0M Q . On the low-mass side (i.e., 
M+ < 2.0-Mq ), the differences are bigger; for ex- 
ample, equation (JTJ) predicts that a cr i t is bigger by 
about 25 % (0.4AU) than its numerical value for 
AI± = 1.5M Q . This is because the stellar mass loss 
is relatively effective in pushing planets outwards 
in the low-mass cases. 

Next, as seen in Fig. [5J there is a sharp transi- 
tion in between 1.7 Mq and 2.1 M Q : a cr it changes 
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Fig. 3. — The critical initial semimajor axis, a cr it, 
inside which planets are engulfed by host stars be- 
fore reaching the end of the helium-burning phase for 
five choices of stellar metallicity, Z±. Dotted, dashed, 
solid, dot-dashed, and double-dot-dashed lines repre- 
sent the results for =0.03, 0.02, 0.01, 0.001, and 
0.0001, respectively. In these calculations, the planet's 
mass is 1 Mj. 



by approximately one order of magnitude in such a 
narrow mass range. The transition corresponds to 
that of i?tip that is related to the thermal state of 
the stellar helium core. A low-mass star of < 2M Q 
has a degenerate helium core after the exhaus- 
tion of the central hydrogen and undergoes the 
helium flash, which causes significant expansion in 
its RGB phase; in contrast, an intermediate-mass 
star of > 2Mq soon starts the 3a burning at its 
center instead of the helium flash, so that i?ti P is 
relatively small. 

Finally, the dependence of a cr it on M*^ is qual- 
itatively the same regardless of stellar metallicity 
(Z*), but a cr it decreases with decreasing Z*, as 
illustrated in Fig. [3] This is also related to the 
dependence of i? t i P on Z*. Because of low opac- 
ity in the atmosphere and of low efficiency of the 
CNO cycle at the hydrogen-burning shell, metal- 
poor RGB stars expand less compared to metal- 
rich RGB stars. 

4. COMPARISON WITH OBSERVA- 
TION 

The critical semimajor axis, a cr it, that we have 
calculated numerically is compared with measured 



Fig. 4. — Comparison between the theoretical sur- 
vival limits and the semimajor axes of known planets 
around clump giants. The solid lines represent the 
critical initial semimajor axes, Ocritj f° r stellar metal- 
licity Z* = 0.01 and planetary mass M v = 1 and 20 
Mj . The semimajor axes of the planets are represented 
by symbols with error bars: filled squares and circles 
correspond to values from the literature and those re- 
evaluated with MESA, respectively. The dashed lines 
represent radial-velocity semi-amplitude of 40 ms -1 
for M p = 1, 2, and 3 Mj. 



semimajor axes of planets so far detected around 
GK clump giants in Figure S) a cr it for M p = 1 Mj 
and 20 Mj and Z* = 0.01 are represented by solid 
lines. The measured semimajor axes are repre- 
sented by symbols with error bars — filled squares 
are data from the literature, while filled circles are 
values that we have re-evaluated using theoreti- 
cal evolution curves generated with MESA. Pub- 
lished values of the stellar masses were eval uated 

20001: 



based on other stellar codes ( Girardi et al 



Lejeune fc Schaererll200ll ) 



Re-evaluation of stellar masses has been done 
just for self-consistency; that never means our val- 
ues are more reasonable than the published ones. 
Because stars of interest in this study are located 
in a region crowded with evolutionary curves in 
the H-R diagram, it is natural that the estima- 
tion of stellar mass includes large uncertainty, 
which is beyond the scope of this study. Note 
that the thre e diffe rent cod es, namely, MESA 
Girardi et al.l ( 2000h 's. and Leieune fc Schaerer 
(1200 lh 's yield similar values of a cr it! as found in 



section E21 The evaluated values of stellar masses 
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are listed on Table [TJ As a result of stellar-mass 
re-evaluation, planetary masses and semimajor 
axes are also changed. These values are also listed 
on Table 1. In comparing observational data with 
theoretical curves, we have assumed that stars are 
located in HeB or post-HeB phases. 

Figure @] indicates two important things: First, 
one finds that almost all the symbols are above 
the curves of a cr it- Although the mean values for 
a few planets are smaller than et cr it, those error 
bars extend out to the curves of a cr it- Second, 
however, one realizes that on the high-mass side 
(i.e., M* > 2.1M Q ), planets exist far above a cr i t . 
There seems to be a gap between the theoretical 
survival- limit and the observed existence-limit of 
planets. Thus, our comparison with observation 
suggests that planet engulfment by host stars may 
not be the main reason for the lack of SPGPs, at 
least, for M+ > 2.1M@. Note that we have to keep 
in mind that the number of planets may be still 
statistically insufficient. 

The presence of the gap is not due to observa- 
tional biases. Clump giants typically show intrin- 
sic radial- velocity variability of a = 10-20 m s _1 
(Sato et al. 2005). Given that we can only detect 
planets that impart radial- velocity semiamplitudc 
of K > 40 m s™ 1 (i.e., 2a) to their host stars, we 
indicate the detection limits by dashed lines for 
M p = 1, 2, and 3 Mj in the figure. While the de- 
tectable limit for planets of < 1 Mj is below a cr it, 
that for planets of > 2 Mj is significantly above 
a„it. Because the masses of the planets shown in 
the figure are all more than 2 Mj , the gap that we 
have found is not due to such observational biases. 

5. SENSITIVITY TO UNCERTAINTIES 



5.1. Tidal dissipation 

For tidal dissipation, this study is based on the 
widespread idea of turbulent viscosity (eqs. [2]- 
[4]). It remains uncertain how to deal with fast 
tide when the planetary orbital period, P, is 
shorter than the eddy turnover timescale, T. For 
such a situation to be taken into account, the fac- 
tor / is added in equation (fTJ) to weaken the stel- 
lar tide. We have adopted the Kolmogorov scal- 
ing, namely, / oc (P/T) 2 , which supposes that 
only small eddies whose intrinsic turnover times 
are shorter than the orbital period contribute 




1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3 3.1 3.2 : 
Host Star's Initial Mass [solar mass] 

Fig. 5. — The impacts of uncertainties in turbulent 
viscosity and eccentricity on the critical initial semi- 
major axes. The three solid lines show the results with 
different frequency-dependent turbulent viscosities, 
namely, the power index, n, of / = min[l, (P/2T)"] 
being 0, 1, and 2. The dashed line corresponds to the 
case with constant quality factor = 10 5 (see text 
for more detail). The dash-dotted line corresponds to 
the case with finite eccentricity of 0.25 (and n = 2). 
In those calculations, = 0.02 and M p = 20Mj. The 
symbols with error bars represent the semimajor axes 
of the planets detected around GK clump giants. 



to th e tidal dissipation (jGoldreich fc Nicholson 
19771 ). However, the Kolmogorov picture is known 
to be too simplified for the turbulent viscosity. 
The dependence of the viscosity on P/T is a 
matt er of debate (e.g., IZahnl 20081 Penev et al 



120071 and references therein) . For example, recent 
hydrodynamic simulations by iPenev et al. ( 2007 ) 
demonstrate that the dependence is close to linear 
(i.e., / oc P/T), which happens to be consistent 
with that proposed by IZahnl ( 19661 ) who assumed 



that the largest eddy dominates tidal dissipation. 

Figure [5] shows a cr it calculated for three differ- 
ent values of the power index n = 0, 1 and 2 (solid 
lines) for / = min[l, (P/T)"}. It also plots the 
measured semimajor axes of the planets of inter- 
est. In calculating o cr it, we assumed the planet's 
mass is 20 Mj, approximately the largest value of 
the minimum masses of planets presented in Fig. [5] 

As n increases, the tidal dissipation becomes 
less efficient, when P < T. In other words, de- 
creasing n moves the curve of a cr it upwards in 
Fig. [SJ As seen in the figure, a cr it for n = 1 is 
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larger by about 50% than that for n — 2 on the 
high-mass side. Note that a cr it hardly differs on 
the low-mass side, because a cr it is so large that P 
is almost always larger than T, namely, / = 1. 
In an extreme case of n = (i.e., / = 1), et cr i t is 
about twice as large as that for n = 2 on the high- 
mass side. Even in this case, decreasing n does 
not fill the gap. 

Finally we comment on the Q-prescription for 
stellar tide. Equation (JTJ can be written in the 
form by the use of the quality factor, Q'+ (e.g., 
Goldreich fc Soterlll966r i. 



1 da 

a dt 



G_ 

Ah 



1/2 



Mp Rl 
OT a 13 / 2 



(8) 



In this study, Q'^ varies with time, while is of- 
ten assumed to be const ant (= 10 5 -10 8 f or MS 
stars in the literature). iNordhaus et al. (2010) 
studied tidal engulfment mainly by AGB stars in 
both cases of variable (/ = 1) and constant 
Q'^ (= 10 5 ). They found that the tidal decay rate 
in the former case is about 4 orders of magnitude 
higher than that in the latter case. In many cases 
of this study, / = (P/2T) 2 <C 1, instead of / = 1. 
Thus, it would be worth checking the outcome for 
constant in our high-mass cases. The dashed 
line in Fig.[5]shows a cr it calculated by equation dHJ 
with Q'^ = 1 X 10 5 . As seen in the figure, the tidal 
decay with Q'+ = 1 x 10 5 yields similar values of 
a„it that we derive with variable Q'+ on the high- 
mass side. This means stellar tide with high 
does not fill the gap neither. 

5.2. Stellar Evolution 

There are several different stellar-evolution 
models in the literature, as mentioned in sec- 
tion [2] We discuss the differences in terms of 
the survival limit, and clarify the cause of the rel- 
ativ ely big difference in q^ it between this study 
and lVillaver fc Livid (|2009l ). 

In Fig. |51 the solid and dashed lines (without 
symbol) represent results that we have obtained 
using MESA with and without overshooting, re- 
spectively. The line with ope n circles is the resu lt 
that we have calculated using ISuda et al.l (|2007l )'s 
code. The two triangles correspond to the val- 
ues of a cr it from Villaver fc Livio ( 20091 ). Those 
two codes do not include overshooting. The lines 
with crosses and squares are the results calcu- 




1.9 2 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 
Host Star's Initial Mass [solar mass] 

Fig. 6. — The critical initial semi-major axes as a 
function of stellar mass calculated with different four 
stellar-evolution co des: The solid and dashed lines are 
from MESA v2258 (|Paxton et al.ll201lh with and with- 
out overshooting, respectively; the lines with crosses, 
with open sq u ares, an d with open circles a re fro m 



iGirardi et al 



( 200Ch' s. iLeieune fc Schaerer! (|200ll )'s 
( 2007lVs codes, respective ly; the filled 
triangles are from IVillaver fc Livid (|2009l ) . The criti- 



and ISuda et al 



cal initial semi-major axes (except the triangles) have 
been calculated until the host star's RGB-tip under 
the assumptions that M cnv = M* and i? olw = in 
the calculations of stellar tide (see eqs [2] and [4]) and 
that Mj, = i n equ ation |T]). As for calculations with 
IGirardi et all (|200d V code, we have assumed that 
is constant with time because values of M* as a func- 
tion of time are unavailable. 



lated with grid data from Girardi et al.l ( 200oT ) and 



Leieune k, Schaererl (|2001l ). respectively, which in 



elude overshooting; those code have been used to 
estimate the mass of planet-host clump giants. 

W hen deriving a rx \t in Fig . [6j except for those 
from Villaver fc Livid ( 20091 ). we have assumed 



that M cnv = M+ and i? C nv = in calculating quan- 
tities relevant to stellar tide, k and T (see eqs. [5] 
and [3]), because no detailed information about 
the stellar internal structure is available. Such 
simplification certainly results in inaccurate eval- 
uation of tidal dissipation (i.e., k/T), but the tidal 
dissipation factor was found to play a minor role 
in a cr it in section 13.21 the most important factor 
is the stellar radius. This means that such simpli- 
fication as to stellar structure does not affect our 
conclusions in this section. 

We have figured out that inclusion of over- 
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shooting probably causes the difference in a cr j t be- 
twe en our calculations (soli d line without symbol) 
and lVillaver fc Livid ( 2009^ 1 's (triangles). Qualita- 
tively, overshooting plays a role in mixing hydro- 
gen and helium, so that stellar models with over- 
shooting contain bigger helium cores compared to 
those without overshooting. This results in higher 
central temperature, which prevents the helium 
core from being degenerated. This means that 
overshooting is effective in moving the sharp tran- 
sition from the high-a cr it to the low-a cr it domains 
leftwards in Fig. |6] That is why we have ob- 
tai ned a value of acrit qu ite different from that 
bv lVillaver fc Livid (j2009h for M* = 2M . 

In any case, the drastic decrease in a cr ;t due 
to avoidance of helium flash does occur, indepen- 
dent of simulation code, although the threshold 
mass somewhat depends on code. Difference in 
code never changes our conclusion concerning the 
survival limit. 

5.3. Eccentricity 

Since we are interested in semimajor axes 
of planets, we have integrated the equation 
for change in semimajor axis (eq. PQ), ignor- 
ing eccentricity, f ollow in g previous studi es (e 



Villaver fc Livid 120091: iNordhaus et all [201 



5: 



However, the planets detected around clump gi- 
ants have finite eccentricities. Here we make a 
simple check of the impact of the finite eccentric- 
ity on the critical semimajor axis. 

The dash-dotted line in Fig. [5] shows a cr it for 
e = 0.25 that we hav e calculat ed by integrating 
the equation (e.g., see Hut 198lh . 



1 da 

a dt 



where 



kMp 
TM+ 



1 



Mp 



/i(e 2 



(J _ e 2)15/2 



(9) 



, / 2% 31 2 255 4 185 R 25 8 , „, 

/i(e 2 ) = + —e 2 + — e 4 + — e 6 + — e 8 . (10) 

The likelihood values of eccentricities of inner 
planets (a < 1.9 AU) are smaller than 0.25. In 
this figure, it is demonstrated that incorporating 
the finite value of e = 0.25 does not affect our con- 
clusion about the presence of the gap. Note that 
the evolution of e itself coupled with that of host 



stars' rotation is an interesting problem, which is 
our future work. 

6. STELLAR-MASS-DEPENDENT PLANET 
FORMATION 

A possibility other than planet engulfment is 
that some stellar-mass-dependent processes hin- 
der the formation of SPGPs around 1.5-3A/© 
stars. For example, ultraviolet measurements 
of young stars indicate that the stellar accretion 
rate (M+) increases rather ra pidly with the stel- 



lar mass, namely, M + oc M£ (jCalvet et al.l 12004 ; 



Muzerolle et al.ll2005l ). This implies that lifetimes 
of protoplanetary disks of heavier s tars a re shorter. 
Based on this fact. lBurkert fc Ida ( 2007 ) supposed 
that protoplanetary disks could dissipate before 
giant planets migrate to the vicinity of the central 
stars in the case of F stars (1.2-1.5Mq). Then, 
they demonstrated that the idea worked to ac- 
count for the observed period valley between the 
hot-Jupiter and cool- Jupiter classe s. Applying the 
idea to heavier stars, Curriei ( 2009f ) proposed that 
SPGPs are rarely formed around ^1.5-3.0 Mq 
stars. 

However, the primordial origin of the paucity of 
SPGPs remains a matter of debate. The negative 
dependence of disk lifetime on stellar mass has not 
yet confirmed theoretically. It is also controversial 
whether planetary migration occurs in disks for 
which such significant stellar accretion is observed. 
Furthermore, there are several competing factors 
that affect the frequency of occurrence of SPGPs 
around ~ 1.5-3 Mq stars. More distant snowlines 
due to high luminosity of such stars also tend to 
reduce the frequency (Ida & Lin 2005). On the 
other hand, more massive disks and shorter Kepler 
periods at given semimajor axes form giant plan- 
ets faster, which results in higher production rate 
of SPGPs and suggests that detection frequency 
of gia nt planets may be higher around high- mass 
stars IjOmiva et al. I l2009h . Further detailed obser- 
vational data on SPGPs will constrain stellar mass 
dependence of planet formation. 

7. SUMMARY AND CONCLUSIONS 

In this study we have derived the survival limit 
beyond which planets are engulfed by their host 
stars in the RGB and HcB phases, by numerically 
simulating both the tidal evolution of planetary or- 
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bits and the evolution of the host stars. Then, we 
have made a detailed comparison of the survival 
limit with observed semimajor axes of planets de- 
tected so far around GK clump giants. The com- 
parison demonstrates that (1) almost all the plan- 
ets are orbiting exterior to our theoretical survival 
limit and (2) the planets around stars with rela- 
tively large masses of > 2.5M Q are far outside the 
survival limit. Those facts suggest that planet en- 
gulfment by host stars may actually happened, but 
it may never be the main reason for the observed 
lack of short-period planets around clump giants. 
Because the number of planets is statistically in- 
sufficient in that high-mass range, we have to say 
that our findings are not definite. The problem 
we addressed have a great influence on our under- 
standing of planet formation. Therefore, further 
surveys for planets around clump giants are highly 
encouraged. 
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Table 1: Properties of stars and planets re- 
evaluated 



Stellar Name M> [Af Q ] a M+ [M ] M p sjnj [Mj] a [AU] Ref. b 



1 1 Pnm c 


o i 4 +0.06 


9 fi +0.4 
Z -°-0.3 


17 1 + 1 - 7 


1 

1 


91 +0.06 
Z1 -0.05 




1 1 TTMi d 
1 ± U 1V11 


-1 o +0.25 
i ' 8 -0.25 


i 7+0.3 


i n s + 1 - 2 

lU.o ^ g 


1 
1 


ci +0.09 
01 -0.09 




14 AnH c 
_L4: rA.llt._l 


2 20 +0 - 13 

z - zu _0 .18 


i 9+0.2 
i - Z -0.3 


o 00+0-4 


u 


fi o+0.03 
08 -0.06 


(A\ (91 


18 Del c 


2 25 +0.06 

z - z -0.05 


2 1+ 01 

-0.1 


9 84 +0 ' 3 

a -°^-0.3 


2 


c-4+0.04 
J -0.04 


(5) (2) 


81 Cet c 


9 .q+0.09 
Z -*°-0.43 


1 6 +0 - 4 

-0.2 


401 +0.6 


2 


17+O.17 
1 ' -0.09 


(4),(2) 


HD 102272 d 


1 q+0-3 

±- -0.3 


1 6 +0 - 8 


5.26+Jj 





t-o+0.08 
00 -0.07 


(6) 


HD 104985 c 


9 i 9 +0.05 
z - lz -0.63 


NA 








(5),(2) 


HD 110014 d 


9 i 7 +0.28 
z ' 1 ' -0.28 


9 -i +0.3 
z ' 1 -0.4 


10.9±i"° 


2 


1 9 +0.09 
iZ -0.15 


(7) 


HD 119445° 


o q+o.4 


q c+0.3 
O.0_ Q 2 


35.0 + ? : ° 


1 


ccc+0.05 
00 -0.03 


(8) 


HD 11977 d 


1 01 +0-21 


1 S+0- 3 

l-S-o.3 


6.29+^ 


1 


on +0.10 
oa -0.11 


(9) 


HD 145457 6 


i q+0-3 

-0.3 


i o +0 - 4 

x- -0.1 


1.89+^ 





fi i +0.08 
-0.02 


(10) 


HD 17092 d 


9 n+0.3 
Z - J -0.3 


1 cc;+0.65 

i.oo _ 025 


3-53lf» 


1 


1 q+0.14 

10 -0.06 


(11) 


HD 173416 c 


2 +0 - 3 


1 O+0.4 

i " 3 -0.4 


2.03 + ° :4 


1 


nn +o.io 

u -Q.ll 


(12) 


HD 180314 6 


2 6 +0 3 

Z -°-0.3 


9 o+O.l 
Z - J -0.4 


20.3 + ™ 


1 


34+0.02 

° -0.08 


(10) 


HD 62509 c 


9 +0.03 
Z ' J1 -0.33 


2 +al 
z - u -o.o 


2 63+° ;i 


1 


fi i +0.03 
01 -0.00 


(13),(2) 


HD 81688 c 


2 07 +0 - 14 

Z ' U ' -0.16 


± - ± -0.2 


1 77 +o.3 

i - ' ' -0.2 





gg+0.05 
uu -0.05 


(5),(2) 


e Tau c 


2 84 +ao7 

Z -°^-0.16 


9 7+O.3 
Z -' -0.2 


7 04+O.5 


1 


qn +0.07 
3 -0.05 


(14),(2) 


£ Aql c 


2 9Q+ - 27 
z - z9 -0.13 


x 4 +0.2 
^ -0.2 


9 Q9+0.2 
Z ' UZ -0.2 





r O+0.02 

00 -0.03 


(5),(2) 


NGC 2423 f 


24 +0.2 
Z '^-0.2 


2 45 + - 35 

Z -^°-0.25 


10 7 +1 ° 


2 


1 1 +0.10 
-0.07 


(15) 


NGC 4349 f 


o n+0.3 

°' y -0.3 


o q+o.5 


19 - 8 -i:? 


2 


qo+0.10 


(15) 



Note. — M*, M p s'mi, and a arc the stellar mass, the planetary minimum mass, and the semimajor axis, respectively. NA 
means that under the assumption that the star is on its HcB or post-HeB, no theoretical track matches its luminosity and 
effective temperature in the H-R diagram. The samples above listed are stars with published values of masses M t > 1.6Mq 
and radii ~ 10 — 20-Rq . When re-evaluating, we have assumed that the solar effective temperature is 5777 K and the solar 
mctallicity is 0.02. 
"Stellar masses, M», in the literature. 

6 (ll ILiu et al.l j200Sf) ; (21 ITakeda et al.l (f200^); (31 iDollinger et afl II2009P) ; (41 ISato et al.l ll2008bH; ( 51 ISato et al.l _2008aH: 
(61 INiedzielski et all feooal; (7) Ide Medeiros et al.| (I2009T ); (81 lOmiva et ail (|200£ ); Ol lSetiawan et al.l d2005l) ; (lOllSato et al" 



2010H ; (111 INiedzielski e t al. [200711; (12l lLiu et all [1200911 ; fl3l lHatzes et allihoog l; (14l lSato et all 1.2007.1 ; (lSl lLovis fe Mavc 



2007) 



c Values of the effective temperature, T e g, and the luminosity, L* and their errors have been taken from the references. 
d Values and their errors of have been derived form those of T e ff and the stellar radius, -R*. 
e Errors of L* have been derived from those of T e ff and R t . 

Stellar masses have been determined so that the stellar ages match the observed cluster age, the stars being assumed to be on 

the HcB. 
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